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Relaxation channels for two-vibron bound states in an anharmonic a-helix protein are studied. 
According to a recently established small polaron model [V. Pouthier, Phys. Rev. E68, 021909 
(2003)], it is pointed out that the relaxation originates in the interaction between the dressed 
anharmonic vibrons and the remaining phonons. This interaction is responsible for the occurrence 
of transitions between two-vibron eigenstates mediated by both phonon absorption and phonon 
emission. At biological temperature, it is shown that the relaxation rate does not significantly 
depends on the nature of the two-vibron state involved in the process. Therefore, the lifetime for 
both bound and free states is of the same order of magnitude and ranges between 0.1 and 1.0 ps 
for realistic parameters. By contrast, the relaxation channels strongly depend on the nature of 
the two-vibron states which is a consequence of the breather-like behavior of the two-vibron bound 
states. 



PACS numbers: 03.65.Ge, 63.20.Ry, 63.22.+m, 87.15.- 
I. INTRODUCTION 

Since the pioneer works of Davydov and co-workers El , 
soliton mechanisms for bioenergy transport in proteins 
have received increasing attention during the last twenty 
five years 0,0 The main idea is that the energy released 
by the hydrolysis of adenosine triphosphate (ATP) can be 
stored in the C=0 vibration (amide-I) of a peptide group. 
The dipole-dipole coupling between peptide groups leads 
to the delocalization of these vibrations and to the forma- 
tion of vibrational excitons, i.e. vibrons. Therefore, the 
strong interaction between the vibrons and the phonons 
of the protein yields a nonlinear dynamics which favors 
the occurrence of the so called Davydov's soliton. 

However, it has been pointed out that the solution of 
the Davydov's problem is rather a small vibron-polaron 
than a vibron-soliton 0, IE IE IE IE El ■ Indeed, the vibron 
bandwidth in proteins is smaller than the phonon cut- 
off frequency so that the non-adiabatic limit is reached. 
During its propagation, a vibron is dressed by a virtual 
cloud of phonons which yields a lattice distortion essen- 
tially located on a single site and which follows instan- 
taneously the vibron (small polaron). Nevertheless, the 
dressing effect leads to an attractive interaction between 
vibrons mediated by virtual phonons. Such an interac- 
tion is responsible for the formation of bound states and 
it has been suggested that proteins can support solitons 
formed by bound states involving a large number of vi- 
brational quanta 0, IE @ • 

Although this formalism gives a comprehensive schema 
for the formation of solitons in proteins, it assumes the 
harmonic approximation for the amide-I vibration. How- 
ever, this approximation failed when several vibrons are 
excited because the intramolecular anharmonicity acts 
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as an additional nonlinear source. As the dressing ef- 
fect, the anharmonicity is responsible for the formation 
of bound states EE El El EE El E3 and the funda- 
mental question of the interplay between both sources of 
nonlinearity has been addressed in a recent paper EE] • In 
this work, we have restricted our attention to the forma- 
tion of two-vibron bound states (TVBS's), only. Indeed, 
although the influence of the anharmonicity in molec- 
ular lattices have been the subject of intense research 
during the last decade, this research was essentially re- 
stricted to classical lattices (for a recent review see for 
instance Refs. [ItI ITil Il9l |). In particular, the formation 
of discrete breathers, i.e. highly localized nonlinear vi- 
brations, has been demonstrated. However, in spite of 
the great interest that these classical nonlinear objects 
have attracted, no clear evidence has been found for their 
existence in real lattices. By contrast, TVBS are have 
been observed in several low-dimensional molecular lat- 
tices [^EimmSHEEIillillil. These quantum 
objects correspond to the first quantum states which ex- 
perience the nonlinearity and can thus be viewed as the 
quantum counterpart of breathers or soliton excitations 
|12| . Their characterization is thus essential and appears 
as a first step to understand the formation of multi- vibron 
solitons. 

In Ref-EH, we have shown that the anharmonicity 
modifies the vibron-phonon interaction which results in 
an enhancement of the dressing effect. Anharmonic vi- 
brons are thus more sensitive to the dressing than har- 
monic vibrons. Moreover, both nonlinear sources favor 
the occurrence of two kinds of bound states which the 
properties strongly depend on the anharmonicity. In 
the harmonic situation, the two bound states appear as 
combinations of states involving the trapping of the two 
vibrons onto the same amide-I mode and onto nearest 
neighbor amide-I modes. By contrast, the intramolecular 
anharmonicity reduces the hybridization between these 
two kinds of trapping so that the low frequency bound 
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state refers to the trapping of the two vibrons onto the 
same amide-I mode whereas the high frequency bound 
state characterizes the trapping onto nearest neighbor 
amidc-I vibrations. 

In this study, the dynamical coupling between the 
dressed anharmonic vibrons and the remaining phonons 
was disregarded. Therefore, the present paper is devoted 
to the characterization of this coupling and to a detailed 
analysis of the relaxation pathways. The TVBS lifetime 
is determined with a special emphasis on the influence of 
the different nonlinear sources. 

At biological temperature, the lifetime of the Davy- 
dov's soliton is still an open question. It has been shown 
that the amide-I excitation, in vivo, corresponds to a lo- 
calized state [2^, • Instead of traveling in a coherent 
manner, it follows a stochastic, diffusional path along the 
lattice. In other words, the single-vibron Davydov soliton 
does not last long enough to be useful at biological tem- 
peratures and it has been shown that two-vibron solitons 
are more stable and appear as good candidates for bioen- 
ergy transport j3J,|33|. However, recent calculations per- 
formed by Ivic et al. clearly show that the multi- 
vibron soliton lifetime is of about a few picoseconds, i.e. 
the same order of magnitude as the single-vibron soliton 
lifetime found by Cottingham and Schweitzer [3^, . 

The paper is organized as follows. In Sec. II, the 
dressed anharmonic vibron point of view described in de- 
tails in Ref. 0] is first summarized. Then, the coupling 
Hamiltonian between these anharmonic polarons and the 
remaining phonons is determined. The TVBS relaxation 
rate is expressed in Sec. Ill and studied numerically in 
Sec. IV. The results are finally discussed in Sec. V. 



II. VIBRON-PHONON HAMILTONIANS AND 
TWO-VIBRON EIGENSTATES 

A. The general vibron-phonon Hamiltonian 

According to the original Davydov's model, the collec- 
tive dynamics of the amide-I vibrations is described by a 
one-dimensional lattice with N sites containing the C=0 
vibrations. The nth amide-I mode is assumed to behave 
as a high frequency anharmonic oscillator described by 
the standard operators 6+ and b n . This oscillator is char- 
acterized by its harmonic frequency ojq and by the cubic 
and quartic anharmonic parameters 73 and 74, respec- 
tively. Finally, the dipole-dipole coupling between near- 
est neighbor amide-I modes is described by the constant 
J. These C=0 vibrations interact with the phonons of 
the lattice which characterize the dynamics of the exter- 
nal motions of the peptide groups. Within the harmonic 
approximation, the phonons correspond to a set of N low 
frequency acoustic modes labeled {q} and described by 
the phonon operators a+ and a q . The frequency of the 
gth mode is defined as Cl q = O c | sin(q/2) |, where fl c 
denotes the phonon cutoff frequency. 

As shown in Ref. |l6j . a unitary transformation is 



performed to remove the intramolecular anharmonicity 
of each amide-I mode and a modified Lang-Firsov trans- 
formation is applied to renormalizc the vibron-phonon 
interaction. As a result, the vibron-phonon Hamiltonian 
is defined as (within the convention K = 1) 



H 



Bb+ +1 b+b n+l b n 



Ji[e+(N n - i)e n+1 (N n+1 )b+b n+1 + h.c] 
J 2 [e+ 2 (v„ - ^)e 2 n+1 (N n+1 + \)b+Hl +1 + h.c] 

J 3 [6+(V„ - l)e n+1 (N n+1 )b+[N n + N n+1 ]b n+1 + H.c] 

J2 n i a t a i CO 



where N n = b+b n and A = 3O7I/W0-674. In Eq.JTJ), H.c. 
stands for the Hermitian conjugate and the different pa- 
rameters are expressed in terms of both the anharmonic 
parameters and the small polaron binding energy Eb as 

d> Q = u -2A-B-(l + 4r])E B 

A = A+{1 + 8ti)Eb:B = B + (1 + 4t 1 )E b 

B = 144J(^) 2 ; Ji = J(l + 44(^) 2 -12^) 



U! 



U) 



J2 = 4J(^) 2 ;J 3 = J(22(^) 2 - 12^i) 

v = 120(— ) 2 - 12— (2) 

Note that the small polaron binding energy Eb and 
the anharmonic parameter A appear as the relevant pa- 
rameters to characterize the nonlinearity of the system. 
In Eq.JU, @ n (N n ) stands for the dressing operator ex- 
pressed as 



B n (N n ) = exp(-Q n [l + 2 ?? + 2 v N n }) 
where Q n is defined as 



(3) 



?n = E 



E B 



sin(q) 



2Nfl q i I sin(q/2) | 



%qn a+ - H.c (4) 



The Hamiltonian Eq. describes the dynamics of an- 
harmonic vibrons dressed by virtual phonons, i.e. an- 
harmonic small polarons. It takes into account on the 
intramolecular anharmonicity up to the second order 
and allows for a renormalization of the main part of 
the vibron-phonon coupling. However, this coupling re- 
mains through the dressing operators O n (V n ) which de- 
pend on the phonon coordinates in a highly nonlinear 
way. Therefore, to separate the vibron degrees of free- 
dom from the phonon coordinates, a mean field proce- 
dure is applied. The full Hamiltonian H is thus writ- 
ten as H = H e ff + Hp + AH, where H p is the phonon 
Hamiltonian and where H e ff =< H > — H p denotes 
the effective Hamiltonian of the dressed anharmonic vi- 
brons. AH = H— < H > stands for the remaining part 
of the vibron-phonon interaction. The symbol < ... > 
represents a thermal average over the phonon degrees of 
freedom at temperature T. 
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B. The effective vibron Hamiltonian and the 
two-vibron eigenstates 

The effective dressed anharmonic vibron Hamiltonian 
is written as 



b+b n - Ab+ 2 bl - Bb+ +1 b+b n+1 b n 



-Ji[#(JV„ + N n+1 )b+b n+1 + H.c] 
-J 2 [$(N„ + N n+1 ) 4 b+ 2 b 2 n+1 + H.c] 
-J 3 [$(N n + N n+1 )b+[N n + N n+1 ]b n+1 + H.c 



(5) 



where $(X) = exp(-S(T)[l + 2rj + 2rjX]) and where 
S(T) is the coupling constant introduced by Ivic and co- 
workers as (fcs denotes the Boltzmann constant) 



S(T) = 



(6) 



In Ref. 0|, a detailed analysis of the two-vibron 
energy spectrum of the Hamiltonian H e ff is p resented. 
Within the number state method dllS El E3, the 
two-vibron wave function is first expanded as | ^) = 
Y^, "2) I Tii, n 2) where | ri\, 712) denotes a local basis 

vector characterizing two vibrons located onto the sites 
ni and n 2 , respectively. Note that the restriction n 2 > rii 
is applied due to the indistinguishable character of the vi- 
brons so that the dimension of the two-vibron subspace is 
iV(iV + 1)/2. Then, by taking advantage of the lattice pe- 
riodicity, the wave function is expanded as a Bloch wave 
as 

*(n 1 ,n 2 =n 1 +m) = -L £ e lk ^ +m ^t> k {m) (7) 



where the total momentum k is associated to the mo- 
tion of the center of mass of the two vibrons whereas 
the resulting wave function ^(771) refers to the degree 
of freedom m which characterizes the distance between 
the two vibrons. Since k is a good quantum number, 
the Hamiltonian H e f f appears as block diagonal and the 
Schrodinger equation can be solved for each k value. For 
a given k value, the protein exhibits (N+l) /2 eigenstates 
I ^fco-), where the index a = 1, (N + l)/2 refers to the 
position of the states along the energy axis. Due to the 
nonlinear sources, there are two different eigenstates, i.e. 
the two-vibron free states (TVFS) and the TVBS. The 
TVFS correspond to a derealization of the separating 
distance m between the two vibrons. The wave function 
\[ / fc(m) behaves as a plane wave and the TVFS belong to 
an energy continuum. By contrast, TVBS correspond to 
a localization of the separating distance between the two 
vibrons and characterize the trapping of the two quanta 
over only a few neighboring sites. We have shown that 
the protein supports two kinds of bound states, called 
TVBS-I and TVBS-II, respectively. The TVBS-I, de- 
noted I ^k,a=i), are located below the TVFS continuum 
over the entire Brillouin zone whereas for TVBS-II, two 



situations occur depending on the strength of the small 
polaron binding energy. For small values of Eb, the band 
disappears inside the continuum when | k | is lower than 
a critical wave vector k c whereas, for strong values of Eb, 
the band is located below the continuum over the entire 
Brillouin zone. As a result, the notation | ^>k.a=2) refers 
either to a free state or to TVBS-II, depending on the 
situation. In the harmonic situation, both TVBS-I and 
TVBS-II appear as combinations of states involving the 
trapping of the two vibrons onto the same amide-I mode 
and onto nearest neighbor amide-I modes. By contrast, 
the intramolecular anharmonicity reduces the hybridiza- 
tion between these two kinds of trapping so that TVBS-I 
refers to the trapping of the two vibrons onto the same 
amide-I mode whereas TVBS-II characterizes the trap- 
ping onto nearest neighbor amide-I vibrations. 



C. The vibron-phonon coupling Hamiltonian 

By comparing Eq.Q arid Eq.JSJl, it is straightforward 
to show that the coupling Hamiltonian AH corresponds 
to a modulation of the different lateral contributions de- 
scribing vibron hops, i.e. the terms proportional to J\, 
J 2 and J3 in Eq.Q). However, in a-helix proteins, it has 
been shown that J 2 ~ J3 ~ Ji/ujq 01- ^ s a result J 2 
and J3 are of about three orders of magnitude lesser than 
Ji and can be neglected. The coupling Hamiltonian AH 
is thus written as 

AH = —Jx [&t(N n + l)e n+s (N n+s )- 

n,S=±l 

(Q+(N n + l)Q n+s (N n+5 ))]b+b n+s (8) 

In addition, the small polaron binding energy is of about 
one order of magnitude smaller than the phonon cutoff 
frequency so that the dressing operator Eq. © can be lin- 
earized H@. As a consequence, by neglecting the rather 
small parameter rj in Eq.Q |16|. the coupling between 
the anharmonic polarons and the remaining phonons is 
finally expressed as 



AH^- &J{n,n + 8)b+b n+s 



n,<5=±l 



where 




(9) 



(10) 



sin(q/2) 



H.c. 



Therefore, within the anharmonic polaron point of view, 
the main contribution of the coupling with the remain- 
ing phonons corresponds to a random modulation of the 
single- vibron hopping constant. As shown in the follow- 
ing sections, this coupling is responsible for the relaxation 
of the two-vibron eigenstates. 
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III. TWO-VIBRON BOUND STATE 
RELAXATION RATE 

Due to the coupling Hamiltonian AH (Eq.@), TVBS 
do not represent exact eigenstates of the whole polaron- 
phonon system. More precisely, this coupling is responsi- 
ble for the occurrence of transitions between two-vibron 
states mediated by the emission or the absorption of 
acoustic phonons. Therefore, by using the Golden rule 
formula, the rate for the transition from a TVBS | ^ka b ) 
with frequency 0Jka b to another state | ^k'a 1 ) with fre- 
quency u>k'<j' is expressed as 

S(uJkcr b + n a — OJk'a' - tya) (11) 

Eq. describes a transition in the course of which the 
phonon bath evolves from an initial state | a) with fre- 
quency f2 Q to a final state | (3) with frequency Qp. Since 
the bath is assumed to be in thermal equilibrium at tem- 
perature T, a statistics is realized over the initial state 
with the probability occupation P a and a sum over all 
the possible final bath states is performed. By inserting 
the expression of the coupling Hamiltonian AH (Eq.JJJJ), 
the total rate for leaving the state | ^k<y b ) obtained by 
summing over all possible transitions is expressed as 

poo 

n,5 n' ,S f k' ,a r 

(AJ(n,n + S,t)AJ{n',n' + 6',0)) (12) 

where (...) stands for an average over the phonon bath 
and where the operators AJ depend on time t accord- 
ing to an Heisenberg representation with respect to the 
phonon Hamiltonian H p . 

As shown in Ea. ll2|) . the TVBS relaxation rate is ex- 
pressed in terms of the Fourier transform of the corre- 
lation function of the coupling AJ. The characteristic 
time of this rate is the correlation time r c of the phonon 
bath which corresponds to the time for which the corre- 
lation functions vanish. In a general way, r c is about 1 
ps for phonons in low-dimensional molecular lattices |l4j . 
We thus assume that this correlation time is sufficiently 
small in order to neglect the spatial correlations in the 
phonon bath. As a result, the correlation functions of 
the coupling A J which appear in Ea. i|12|) are nonzero if 
n = n'+8' and n' = n+S, only. Therefore, by performing 
both the time integration as well as the thermal average 
in Ea. l|12|l . the relaxation rate is finally expressed as 



W ka 



[F(uk'a> ~ Ukcr b )n(uJk>a> - WfecJ + 

F(u kab - Wfev'X 1 + n(uj kcrb - uJk'a'))} (13) 



where n(O) denotes the Bose-Einstein phonon distribu- 
tion at temperature T and where the coupling distribu- 
tion F(Q,), which measures the probability for the ex- 
change of a phonon with frequency Q during the process, 
is defined as 



F(n) 







(W-J 2 if ^ > 

if n < o 



(14) 



In Ea. (|13fl . Zka^k'cr' characterizes the strength of the 
coupling between the two-vibron eigenstates | ^ka) and 
I ^k'a') due to the vibron-phonon interaction. This cou- 
pling is expressed as 



Z, 



ka — >k' 



n,5=±l 



(* feCT | b+b n+s | * feV ') I 2 (15) 



After some algebraic manipulations, Zka^k'a' is ex- 
pressed in terms of the wave functions as 



->ko— >k'a' 



1_ | \ - ^ka{m) ^k'a'jm- 1) i{ k'-k) m /2 

N 1 ^ A(m) [ A(m- 1) 

2 



A(m + 1) 1 

1 . V- ^k'a'jm) fr fc<T (m- 1) l(k '-k)r. 

N 1 ^ A(m) 1 A(m- 1) 

^ka{m+ 1) l(k '-k),n/2 |2 

A(m + 1) 1 



;/2 



(16) 



where the convention ^ka(~ 1) is used. 

As shown in Ea. (|13fl . the rate depends on the temper- 
ature through the average number of phonons. More- 
over, the temperature is involved in the definition of the 
two-vibron wave functions due to the dressing effect |16| . 
The relaxation rate exhibits two contributions connected 
to the absorption (term proportional to n(fl)) and to the 
emission (term proportional to 1 + n(Q)) of an acoustic 
phonon, respectively. Note that Eq. l|13H clearly shows 
that the relaxation rate Wk ab is expressed as the sum 
over the rate connected to the different relaxation chan- 
nels Wka b ^k' a' ■ As a result, each channel can be charac- 
terized separately. 

At this step, the diagonalization of the Hamiltonian 
H e fj realized in Ref.fla] allows us to calculate both 
the two-vibron eigenstates | ^ka-} and eigenenergies oJka- 
Then, by using Eqs.lplll-lpjjl, tnc TVBS relaxation rate 
can be computed. This procedure is illustrated in the 
following section. 



IV. NUMERICAL RESULTS 

In this section, the previous formalism is applied to 
compute the TVBS relaxation rate in an anharmonic 
a-helix protein. The intramolecular anharmonicity is 
described by a single parameter, namely the anhar- 
monic constant A, which ranges between and 10 cm -1 
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[l6l l35l l3^ |. The small polaron binding energy Eb is 
taken as a parameter which extends from to 15 cm . 
The phonon cutoff frequency f2 c is fixed to 100 cm" 1 and 
the hopping constant is set to J = 8 cm -1 . 

The temperature dependence of the zero wave vector 
TVBS-I relaxation rate (full circles) is shown in Figs, 
la and lb for two typical situations. The empty circles 
correspond to the rate for the relaxation over all the other 
TVBS-I whereas the empty squares represent the rate for 
the decay into the set of the second eigenstates | ^ka=2)- 
Note that as remained in Sec. II. B., such states refer to 
either free or bound states (TVBS-II) , depending on the 
nonlinearity. 

When A=8 cm -1 and E B — 4 cm -1 (Fig. la), the 
relaxation rate exhibits a quasi-linear dependence ver- 
sus the temperature, excepted at low temperature. More 
precisely, the linear regime is reached when the tempera- 
ture is greater than 50 K whereas the rate shows a power 
law dependence at low temperature. The relaxation rate 
is equal to 0.049 cm -1 at T — 5 K and reaches 7.64 
cm -1 at T = 315 K. At low temperature, Fig. la clearly 
indicates that the main mechanism for the relaxation in- 
volves the decay of the zero wave vector TVBS-I into 
the other TVBS-I. For instance, this channel represents 
99.8 % of the relaxation at T — 5 K. As increasing the 
temperature, the relaxation over the other TVBS-I de- 
creases and the rate for the relaxation over the states 
I ^k<T=i) increases very slightly. Indeed, at T = 315 K, 
the relaxation over the other TVBS-I represents 55.54 %. 
However, the second channel, i.e. the relaxation over all 
the states a = 2, represents only 7.22 % which indicates 
that 37.24 % of the relaxation involves the decay of the 
TVBS-I into the TVFS continuum. 

When Eb = 12 cm -1 (Fig. lb), the relaxation rate for 
the zero wave vector TVBS-I exhibits almost the same 
temperature dependence as in the previous case. Never- 
theless, the rate is more important since it is equal to 0.19 
cm" 1 at T = 5 K and reaches 19.58 cm" 1 at T = 315 K. 
However, the main difference with the previous case orig- 
inates in the nature of the relaxation channels. Indeed, 
although the relaxation over the other TVBS-I remains 
the main pathway at low temperature, this is no longer 
true at high temperature. Indeed, as shown in Fig. lb, 
the decay into the other TVBS-I represents almost 100 
% of the relaxation at T = 5 K. However, as increasing 
the temperature, the relaxation over the other TVBS- 
I strongly decreases whereas the rate for the relaxation 
over the states | ^!ka=2) increases and becomes the dom- 
inate contribution (the transition occurs around 130 K). 
At high temperature, i.e. T — 315 K, the relaxation ac- 
cording to the second channel represents 85.60 % whereas 
the decay into the other TVBS-I charaterizes 10.84 % of 
the global rate. As a consequence, for such a strong non- 
linear situation, the decay of the TVBS-I into the TVFS 
is no more than 5 % at high temperature. 

The behavior of the TVBS-I relaxation rate as a func- 
tion of the anharmonicity is displayed in Figs. 2. The cal- 
culations are performed at T — 310 K and for three dif- 



ferent values of the small polaron binding energy. When 
Eb = 4 cm -1 (Fig. 2a), the rate slightly decreases as 
the anharmonicity increases. It is equal to 8.28 cm -1 
when A = and to 7.26 cm -1 when A = 10 cm -1 . In 
marked contrast, the rate for the decay into the other 
TVBS-I first increases to reach a maximum equal 6.02 
cm -1 when A = 2 cm -1 . Then, it decreases and is equal 
to 3.48 cm -1 when A=10 cm -1 . As shown in Fig. 2a, the 
rate for the decay into the second eigenstates | ^ka=2) 
is rather small whatever the anharmonicity although it 
increases when A increases. Consequently, in this low 
nonlinear regime {Eb — 4 cm -1 ), the main part of the 
relaxation of TVBS-I involves the decay into both the 
other TVBS-I and the TVFS continuum. 

As when increasing the small polaron binding energy, 
(Fig. 2b and 2c), the TVBS-I relaxation rate behaves 
in a similar way with respect to the anharmonicity and 
slightly decreases as A increases. However, the rate in- 
creases with Eb since it is equal to 16.31 cm -1 when 
E B = 8 cm -1 and A = (Fig. 2b) and reaches 22.73 
cm -1 when Eb — 12 cm -1 and A — (Fig. 2c). In 
a marked contrast, the relaxation pathways are strongly 
modified when the small polaron binding energy is in- 
creased. Indeed, when E B = 8 cm -1 (Fig. 2b), the 
rate for the decay into the other TVBS-I decreases as A 
increases. By contrast, the rate for the decay into the 
second eigenstates | ^ka=2) increases. This second chan- 
nel becomes slightly more efficient than the first channel 
when the anharmonicity exceeds 7 cm . When A = 10 
cm" 1 , the first channel represents 21.92 % of the relax- 
ation whereas the second channel characterizes 43.31 % 
of the decay. Such a behavior appears more pronounced 
when the small binding energy is set to Eb — 12 cm" 1 . 
In that case, Fig. 2c clearly shows that the decay of the 
TVBS-I into the TVFS continuum is rather weak. Its 
contribution is less than 5 % when A is greater than 3.5 
cm" 1 . For a small anharmonicity, the rates for the decay 
into the other TVBS-I and into the second eigenstates 
are of the same order of magnitude. However, as in- 
creasing the anharmonicity, the second channel becomes 
dominant since the corresponding rate represents almost 
90 % of the global rate when A=10 cm" 1 . 

The dependence of the TVBS-I relaxation rate on the 
small polaron binding energy is shown in Fig. 3 for 
T = 310 K, and A = 8 cm" 1 . The global rate evolves 
in a quasi-linear way and varies from 1.98 cm -1 when 
E B = 1 cm" 1 to 23.19 cm" 1 when E B = 15 cm" 1 . As 
shown in Fig. 3, the more surprising results correspond 
to the behavior of the rates connected to the first and 
to the second channel. Indeed, for small Eb values, the 
decay into the other TVBS-I is the dominant relaxation 
pathway. For instance, when Eb — 2 cm" 1 , this first 
channel represents 67.66 % of the relaxation whereas the 
contribution of the second channel is 2.13 %. Therefore, 
30.21 % of the relaxation involve the decay of the TVBS-I 
into the TVFS continuum. However, as Eb increases, the 
rate connected to the second channel increases and be- 
comes the main contribution for strong Eb values. When 
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Eb = 14 cm -1 , the second channel represents 91.54 % of 
the relaxation whereas the contribution of the first chan- 
nel is 6.81 %. Note that both channels contribute in a 
similar way when Eb is about 7.5 cm -1 . 

The correlation between the relaxation channels and 
the nature of the two-vibron eigenstates, is illustrated in 
Figs. 4 for T = 310 K and A = 8 cm -1 . The upper panel 
represents the corresponding two-vibron energy spectrum 
whereas the lower panel displays the wave vector depen- 
dence of the relaxation rates. More precisely, open circles 
characterize the rate Wo,i->ko-=i f° r the decay of the zero 
wave vector TVBS-I into the TVBS-I with wave vector k. 
In the same way, open squares correspond to the decay 
of the zero wave vector TVBS-I into the state | ^ka=2)i 
i.e. Wo,i->ko-=2- 

When E B = 4 cm- 1 (Fig. 4a), the TVBS-II band 
occurs at the end of the first Brillouin zone, only. For 
the first relaxation channel, i.e. the relaxation over the 
other TVBS-I, the rate decreases as the modulus of the 
wave vector increases. In other words, the decay into low 
wave vector TVBS-I is the dominant relaxation pathway. 
Although such an effect is not correlated to the energy 
spectrum, this is no longer true for the k dependence of 
the rate connected to the second channel. Indeed, Fig. 
4a clearly shows that the second channel opens at the end 
of the first Brillouin zone, only, i.e. when the TVBS-II 
band occurs. Therefore, this second channel becomes the 
dominant relaxation pathway at the end of the Brillouin 
zone. When E B — 10 cm -1 (Fig. 4b), the results are 
slightly different since Eb is strong enough so that the 
TVBS-II band is localized below the continuum over the 
entire Brillouin zone. As a result, the second channel is 
clearly the main mechanism for the relaxation whatever 
the value of the wave vector. 

In Fig. 5, the influence of the small polaron binding 
energy on the relaxation rate of the | ^fc=o,<r=2) eigen- 
state is shown at biological temperature (T=310 K) and 
for A = 8 cm -1 . As in the previous figures, full circles 
characterize the global rate, empty circles represent the 
rate for the decay into the TVBS-I and empty squares 
correspond to the rate for the decay into all the other 
I ^ f fc,<T=2) eigenstates. In addition, empty triangles char- 
acterize the rate for the decay into the TVFS continuum. 
The figure clearly shows that the system exhibits two 
regimes depending on the value of Eb- 

For the small values of Eb, i.e. typically Eb < 8.5 
cm , the rate increases in a quasi-linear way with Eb- 
It is equal to 1.65 cm -1 when Eb = 1 cm -1 and reaches 
13.67 cm -1 when Eb — 8.0 cm -1 . As shown in Fig. 
5, the relaxation into both the TVBS-I and the other 
I ^k,a=2) eigenstates can be neglected. In other words, 
the main relaxation pathway corresponds to the decay 
of the I ^0^=2) eigenstate into the TVFS continuum. 
When Eb becomes greater than a critical value, the 
global rate behaves in a different manner with respect 
to the small polaron binding energy. It increases more 
rapidly than the previous linear regime to reach 38.75 
cm -1 for Eb = 15 cm -1 . Moreover, the rate for the re- 



laxation into the TVBS-I increases with Eb whereas the 
rate for the relaxation into TVFS slightly decreases. For 
strong Eb, the decay into TVBS-I represents 64 % of the 
relaxation whereas the decay into TVFS corresponds to 
33 %. Note that although the rate for the decay into the 
I ^k,a=2) eigenstates sligtly increases around the transi- 
tion, it decreases as Eb increases and can be neglected 
for a strong nonlinearity. 

Finally, the relaxation rate of the | ^k=Q,a=2) eigen- 
state versus the anharmonicity is shown in Fig. 6 for 
T=310 K and for E B = 10 cm" 1 . As in Fig. 5, the 
different rates exhibit two regimes depending on the an- 
harmonicity. For a small anharmonicity, i.e. A < 4 cm , 
the global rate is almost independent on the anharmonic 
parameter and it is equal to 17.00 cm -1 . In addition, 
the decay into the TVFS continuum appears as the main 
mechanism for the relaxation. In a marked contrast, for a 
strong anharmonicity, i.e. A > 4 cm -1 , the rate increases 
as the anharmonic parameter increases and reaches 23.19 
enr 1 when A — 10 cm -1 . The rate for the relaxation 
over the TVFS continuum decreases whereas the rate for 
the decay into the TVBS-I increases. This latter rate 
becomes the most important when A = 10 cm -1 so that 
the decay into TVBS-I represents almost 50 % of the 
relaxation. 



V. DISCUSSION 

To discuss and interpret the previous numerical results, 
let us first consider the behavior of the TVBS relaxation 
rate at low temperature. Since the zero wave vector 
TVBS-I lies at the bottom of the two-vibron energy spec- 
trum, its decay involves phonon absorption, only. The 
TVBS-I relaxation rate is thus proportional to the Bose- 
Einstein distribution which selects the frequency range 
of the phonons which are exchanged. In that context, 
transitions involving low frequency phonons takes place 
at low temperature. Therefore, when the thermal energy 
ksT is lower than the energy gap between bound and 
free states, the TVBS-I can just decays into the other 
TVBS-I, as shown in Fig. 1. Note that phonon emission 
participates in the decay of the other two-vibron states 
so that the corresponding rate reaches a finite value at 
zero temperature (not considered in the numerical anal- 
ysis). As a consequence, the low temperature behavior 
of the rate does not depend on the system nonlinearity 
and essentially originates in the shape of the phonon dis- 
tribution. 

This is no longer the case at biological temperature 
for which an approximate expression of the relaxation 
rate can be determined. To proceed, we assume that the 
two-vibron bandwidth is smaller when compared with 
both the thermal energy ksT and the phonon cutoff 
energy M7 C . Therefore, the Bose-Einstein distribution 
can be linearized according to the temperature and the 
distribution function F(f2) (Eq.JTJJ) can be written as 
F(fl) w f2/fi c . As a result, the relaxation rate Ea. H13|) 
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connected to the two-vibron eigenstate | ^ka) can be ap- 
proximated as 

At this step, Ea. (|17(l can be simplified because the sum 
over k'a' leads to the occurrence of the closure relation. 
Therefore, by using the following identity 

I b+b n b+b n , | V k<7 ) = I | » - n' |) | 2 (18) 

the relaxation rate is finally expressed as 

W k ^ l2 * J ^ kBT {l + \\*Ul)\ 2 ) (19) 

Ea. l|19|) yields a rather good approximation which al- 
lows us to interpret and to understand the numerical re- 
sults. First of all, it accounts for the observed tempera- 
ture dependence of the relaxation rate. Indeed, Ea. l(T^|l 
clearly shows that the rate increases in a linear way as 
the temperature increases, in a perfect agreement with 
the numerical results displayed in Figs. 1. This feature 
originates in the linearization of the Bose-Einstein dis- 
tribution. Note that the wave function ^^(m) depends 
on the temperature in a complicated manner through the 
dressing effect 16] . However, our results indicate that 
such a dependence remains rather small when compared 
with the influence of the Bose-Einstein factor. 

Then, the wave function dependence of the relaxation 
rate allows us to understand the influence of the in- 
tramolecular anharmonicity as shown in Figs. 2 and 
6. For TVBS-I, Figs. 2 clearly show that the rate de- 
creases as the anharmonicity increases. In fact, since 
TVBS-I refer to the trapping of the two vibrons onto 
the same amide-I vibration, the wave function 4 , fc cr= i(m) 
is maximum for m = and decreases with m accord- 
ing to a quasi-exponential way (lfjj. As a consequence, 
when the anharmonicity increases, the trapping process 
is enhanced so that the extension of the wave function 
around m = is reduced. Therefore, | ^ , / C(T= i(0) | 2 in- 
creases whereas | \P fccr=i(l) | 2 , as the relaxation rate, de- 
creases. In Fig. 6, the anharmonicity dependence of the 
rate for the | ^>k=0a=2) eigenstate exhibits two regimes 
which originates in the nature of the state itself. Indeed, 
as pointed out in Sec. II. B., | ^f k=0a=2 ) refers either to 
a free or to a bound state (TVBS-II), depending on the 
nonlinearity. Such a behavior is displayed in Fig. 6 since 
the change of regime corresponds to the transition from a 
TVFS to a TVBS-II. As shown in Ref. [3, | %= 0ct =2) 
corresponds to a TVFS for a small anharmonicity. In 
that case, the probability | ^kcr=2(m) | 2 is almost inde- 
pendent on the anharmonicity and scales as 1/N. Con- 
sequently, the relaxation rate does not depend on A, as 
shown in Fig. 6. By contrast, as when increasing the 
anharmonicity, the state | ^k=o<r=2) becomes bounded. 
It refers to a TVBS-II which characterizes the trapping 



of the two vibrons onto two nearest neighbor amide-I vi- 
brations. Therefore, when the anharmonicity increases, 
the trapping process is enhanced so that the wave func- 
tion , 5fc ( j = 2(?Tt) tends to localize around m = 1. The 
corresponding rate suddenly increases with A as clearly 
shown in Fig. 6. 

Finally, the approximated expression of the rate gives 
a comprehensive explanation of the influence of the small 
polaron binding energy (see Figs. 3 and 5). Indeed, 
Ea. (|19|) clearly shows that the dependence of the rate 
with respect to Eb is twofold. First, the rate depends 
linearly on Eb- This feature originates in the fact that 
the rate is proportional to the intensity of the coupling 
between the dressed anharmonic vibrons and the remain- 
ing phonons (see Ea. ljTTl) '). Then, the Eb dependence is 
included in the wave function dependence of the rate. As 
when increasing the small polaron binding energy, the 
trapping process involved in the formation of the TVBS- 
I is enhanced so that 4 , fc cr= i(l) decreases. Therefore, the 
linear evolution of the TVBS-I relaxation rate with re- 
spect to Eb is slightly damped, i.e. the rate evolves more 
slowly than the corresponding linear law (see Fig. 3). As 
previously, Fig. 5 clearly shows that the nature of the 
state | ^k=0a=2) exhibits a transition as a function of 
the small polaron binding energy. For small Eb values, 
this state is a TVFS which the wave function does not 
significantly depend on the nonlinearity. Therefore, the 
rate evolves in a linearly way versus Eb- By contrast, 
when the small polaron binding energy is sufficiently im- 
portant, this state becomes a TVBS-II |lj|. As a con- 
sequence, the trapping onto two nearest neighbor sites 
is enhanced so that 5 , fe (T= 2(l) increases to reach unity. 
The corresponding relaxation rate suddenly evolves more 
rapidly than the previous linear law (see Fig. 5). 

As shown in Ea. (|19fl . the dependence of the relaxation 
rate on the nature of the corresponding two-vibron state 
arises from the factor (1+ | "$>ka{^) | 2 /2). As a conse- 
quence, the global rate does not depend significantly on 
the fact that the state is a TVFS or a TVBS. In other 
words, the lifetime of both kinds of states is about the 
same order of magnitude. Indeed, | ^fe (T (l) | 2 refers to 
the probability to find the two vibrons onto two near- 
est neighbor sites. In a TVFS, this probability scales as 
1/7V and vanishes asymptotically as TV tends to infinity. 
For a TVBS, such a probability depends on the nature of 
the bound between the two vibrons and ranges between 
and 1. For instance, for a strong nonlinearity, i.e. for 
strong values of both A and Eb, TVBS-I refer to the 
trapping of the two vibrons onto the same amide-I vibra- 
tion whereas TVBS-II characterize the trapping onto two 
nearest neighbor sites. Therefore, | ^^(l) | 2 vanishes for 
TVBS-I whereas it is equal to unity for TVBS-II. In that 
context, if we define W = 128JlE B kBT/nl, the relax- 
ation rates are expressed as Wtvbs-i ~ Wtvfs ~ W 
and Wtvbs-ii ~ 3VF/2. 

In a marked contrast, the nature of the relaxation 
channels drastically depends on the characteristic of the 
two-vibron states involved in the process. Such effects 
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originate in the coupling between the two-vibron eigen- 
states mediated by the vibron-phonon interaction and 
characterized by the constant Z ka ^ k i a i (see Eas. (ll5[) - 
(|16fl % ). As shown in the definition of the coupling Hamil- 
tonian AH (Eq.©), this interaction characterizes tran- 
sitions between two-vibron states in the course of which a 
single vibron is transfered from a given site to its nearest 
neighbor. 

To understand this feature, let us first consider the re- 
laxation channels connected to the decay of the TVBS- 
I. Within the strong nonlinear limit, the TVBS-I wave 
function is essentially localized around m = so that 
v E'fecr=i(m.) is almost equal to 5„ lt Q. Therefore, it is 
straightforward to show that the coupling constant can 
be approximated as Z krT= i^ k '<r' » 4 | fy ka=1 (0) | 2 | 
"ffc'cr' (1) | 2 /N. The corresponding relaxation rate is thus 
written as 

W ka=1 ^ k , a , « ^ | *to=i(0) | 2 | * fe v'(l) I 2 (20) 

The relaxation channel for the TVBS-I corresponds to 
the decay into states which favor the trapping of the two 
vibrons onto two nearest neighbor amide-I vibrations, i.e. 
states for which | 4 r fc' cr '(l) | 2 is maximum. Within the 
strong nonlinear limit, such states refer to the TVBS-II. 
As a consequence, in a perfect agreement with the numer- 
ical results displayed in Figs, la, 2c, 3 and 4b, Ea. ij2U|) 
clearly shows that the decay into TVBS-II is the main 
mechanism for the TVBS-I relaxation. Note that both 
| ^ka=i{0) | 2 and | 4'fc (T= 2(l) | 2 are almost k independent 
so that the rate for the decay over all the TVBS-II states 
is equal to W. However, when the nonlinearity slightly 
decreases, i.e. when either (or both) A and Eb decreases, 
the TVBS-I wave function delocalizes so that ^fc CT =i(0) 
decreases and, in the same time, ^ ka= \(l) increases. As 
a result, the coupling between a given TVBS-I and the 
other TVBS-I is turned on which opens the correspond- 
ing relaxation channel. In the same way, the TVBS-II 
wave function extends itself around m = 1 leading to the 
decrease of vE , fe (T =2(l)- As a consequence, the coupling 
between TVBS-I and TVBS-II decreases, as observed in 
the numerical results. 

Within the strong nonlinear limit, the TVBS-II wave 
function is essentially localized around m = 1 so that 
v I'fcc7=2(™) is almost equal to 5 m ^\. Therefore, by fol- 
lowing the same procedure as in the previous para- 
graph, it is straightforward to show that the coupling 
constant is Z ka=2 ^ k , a , w 4 | * feCT=2 (l) | 2 | * fe ' CT /(0) + 
#fc/ CT /(2)exp(£(fc ~ k'))/V2 | 2 /N. This expression indi- 
cates that TVBS-II decay into both TVBS-I and TVFS, 
but cannot relax over the other TVBS-II. As a result, the 
rate for the decay of a given TVBS-II into a particular 
TVBS-I is approximately given by 

W kG ^ k , a=x « ^ | *to=2(l) | 2 | *fcv=i(0) | 2 (21) 

At this step, by performing the sum over the wave vector 
k' in Ea.l|!n)l. the rate for the decay over all the TVBS-I 



is equal to W. Since the global relaxation rate for the 
TVBS-II is equal to 3 W/2, the rate for the relaxation over 
the TVFS continuum is equal to W/2. These results are 
in a rather good agreement with the numerical results 
shown in Fig. 5. However, when the nonlinearity slightly 
decreases, both the TVBS-I and TVBS-II wave functions 
delocalize so that \E , fe (T= 2(l) and ^/fccr=i(0) decrease. In 
that context, the relaxation rate for the decay of TVBS- 
II into TVBS-I decreases, as shown in Figs. 5 and 6, 
whereas the rate for the relaxation over TVFS increases. 
In the same way, the relaxation channel connected to the 
decay into the other TVBS-II opens because the wave 
function ^ ka=2{fn) does not vanish any more for m = 
or m = 2. 

At this step, let us mention that the previous results 
address the fundamental question of the manifestation of 
the quantum localization, i.e. of the existence of quan- 
tum breathers. Indeed, classical breathers describe time 
periodic and spatially localized nonlinear vibrations. By 
contrast, due to the translational invariance, the center 
of mass of the two vibrons in a TVBS is fully delocalized 
according to a Bloch wave. However, the localized nature 
of a TVBS arises due to the trapping of the two quanta 
around few sites. In that context, although a TVBS does 
not describe a spatially localized field, its degree of lo- 
calization manifests itself through particular correlation 
functions |37J, |38| . Our results clearly establish that such 
correlation functions are involved in the calculation of the 
relaxation rate (see for instance Eq. As a conse- 

quence, the dependence of the relaxation channels versus 
the nature of the two-vibron states can be viewed as a 
manifestation of the quantum localization of TVBS. 

To conclude, let us note that the knowledge of the 
relaxation rates allow us to characterize the population 
dynamics of the two-vibron eigenstates in real time. To 
illustrate this feature, let us consider the strong nonlinear 
limit in which the relaxation schema is rather simple. 
Indeed, the rates are almost wave vector independent so 
that the system can be modeled as a three levels system 
formed by the TVBS-I band, the TVBS-II band and the 
TVFS band (see Fig. 7). The first relaxation channel 
corresponds to the decay of the TVBS-I band into the 
TVBS-II band according to the rate W. The relaxation of 
the TVBS-II band exhibits two main channels connected 
to the decay into the TVBS-I band with the rate W and 
to the decay into the TVFS band with the rate W/2. 
Finally, TVFS are allowed to relax into the other TVFS, 
only. Let Pj , Pjj and Pp denote the populations of these 
bands. Therefore, the time evolution of these populations 
is governed by a master equation expressed as 

Pi = -WPi + WPn 
3W 

Pjj = WPl-—P n 

W 

Pf - -j Pi 1 (22) 

The system Ea. (|22|) . which can be solved straightfor- 
wardly, allows us to follow the relaxation pathways in 
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real time. For instance, if we assume that the system is 
initially set in the TVBS-I band, Ea. (|22|) clearly shows 
that the population Pi(t) decays so that the population 
of the TVBS-II band increases. However, as Pn(t) in- 
creases with time, it tends to relax into the TVFS contin- 
uum. Therefore, we observe the relaxation of the TVBS- 
I band into the TVFS continuum through the TVBS-II 
band. After a time t* — 21.5/W, both Pi and Pjj can 
be neglected whereas Pp is equal to 0.99. At biological 
temperature (T = 310 K), with the realistic parameters 
A = 8 cm -1 , Eb = 12 cm -1 , and J = 8 cm -1 we obtain 
W = 21.36 cm- 1 so that t* = 5.3ps. 

Note that in a recent paper devoted to the TVBS re- 
laxation in a molecular nanowire, it has been pointed out 
that, for a one-dimensional phonon bath, a given TVBS 
essentially decays into the other TVBS |l4|. The origine 
of the different with the present work is twofold. First, 
in Ref.0], the vibron dynamics was described accord- 
ing to a Hubbard model for bosons which yields a single 
TVBS band (almost identical to the TVBS-I). Then, the 
coupling with the thermmal bath was assumed to be re- 
sponsibel for a random modulation of the frequency of 
each molecule, i.e. AH « ^ Auj n b^b n , in marked con- 
trast with the modulation of the vibron hopping constant 
considered in this paper (see Eq.©). 

To summarize, the present paper was devoted to the 
characterization of the two- vibron relaxation mechanisms 
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Figure Caption 

Figure 1 : Temperature dependence of the zero wave 
vector TVBS-I relaxation rate (full circles) for A = 8 
cm -1 and for E B = 4 cm -1 (a) and E B = 12 cm -1 (b). 
Empty circles correspond to the rate for the relaxation 
over all the other TVBS-I whereas empty squares rep- 
resent the rate for the decay into the set of the second 
eigcnstates | , t>ka=2) (see the text). 

Figure 2 : TVBS-I relaxation rate (full circles) versus 
the intramolecular anharmonicity at T — 310 K for for 
E B = 4 cur 1 (a), E B = 8 cm" 1 (b) and E B = 12 cm" 1 
(c). Empty circles correspond to the rate for the relax- 
ation over all the other TVBS-I whereas empty squares 
represent the rate for the decay into the set of the second 
eigenstates | ^ku=2) (see the text). 

Figure 3 : TVBS-I relaxation rate (full circles) versus 
the small polaron binding energy at T — 310 K for for 
A = 8 cm -1 . Empty circles correspond to the rate for 
the relaxation over all the other TVBS-I whereas empty 
squares represent the rate for the decay into the set of 
the second eigenstates | \f r fc<r=2) (see the text). 



Figure 4 : Correlations between the relaxation chan- 
nels and the nature of the two-vibron eigenstates for 
T = 310 K, A = 8 cm" 1 and E B = 4 cm" 1 (a) and 
E B = 10 cm -1 (b). The upper panel represents two- 
vibron energy spectrum whereas the lower panel displays 
the wave vector dependence of the relaxation rates. Open 
circles characterize the rate W / o,i-»/ SCT =i whereas open 
squares correspond to the rate Wo,i-»fc<r=2- 

Figure 5 : Relaxation rate (full circles) of the | 
*k=o,<7=2) eigenstate versus E B at T=310 K and for 
A = 8 cm -1 . Empty circles represent the rate for the 
decay into the TVBS-I and Empty squares correspond 
to the rate for the decay into all the other | ^k,a=2) 
eigenstates. Empty triangles characterize the rate for 
the decay into the TVFS continuum. 

Figure 6 : Relaxation rate (full circles) of the | 
v I'fe=o,(T=2) eigenstate versus A at T=310 K and for E B = 
10 cm' 1 . Empty circles represent the rate for the decay 
into the TVBS-I and Empty squares correspond to the 
rate for the decay into all the other | ^k,<y=2) eigenstates. 
Empty triangles characterize the rate for the decay into 
the TVFS continuum. 

Figure 7 : Three levels model for the population dy- 
namics of the two-vibron states within the strong nonlin- 
ear limit. 
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